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We have performed high-resolution two-photon dark-state spectroscopy of an ultracold gas of 
87 Rb2 molecules in the a 3 Sj state at a magnetic field of about 1000 G. The vibrational ladder 
as well as the hyperfine and low-lying rotational structure is mapped out. Energy shifts in the 
spectrum are observed due to singlet-triplet mixing at binding energies as deep as a few hundred 
GHzx/i. This information together with data from other sources is used to optimize the potentials 
of the a 3 Ej and X 1 E^" states in a coupled-channel model. We find that the hyperfine structure 
depends weakly on the vibrational level. This provides a possible explanation for inaccuracies in 
recent Feshbach resonance calculations. 

PACS numbers: 42.62.FI, 33.20.-t, 37.10.Mn 
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I. INTRODUCTION 



Recently, optical schemes have been developed to se- 
lectively produce cold and dense samples of deeply bound 
molecules in well defined quantum states [1H3 ( see a l so 
d, Q). This has opened up new possibilities for cold 
collision experiments, ultracold chemistry [fHH, and for 
testing fundamental laws via precision spectroscopy 0- 
[Tl| . For such future experiments it is mandatory that the 
location and properties of the available molecular quan- 
tum states are well known and understood. 

Until recently, the a 3 E+ triplet potential of the 
Rb2 molecule has remained largely experimentally un- 
explored. Natural samples of Rb2 are normally found in 
their X 1 !]^ singlet state from which the lowest triplet 
state is somewhat difficult to reach in an optical Raman 
process due to a change in u/g symmetry. Here we report 
precision spectroscopy of 87 Rb2 molecules in this lowest 
triplet state, a 3 S+ (55i/ 2 + 55*1/2), where we resolve 
vibrational, rotational, hyperfine and Zeeman structure 
with an accuracy as high as 30 MHz. We use dark- 
state spectroscopy where a gas of weakly bound Feshbach 
molecules is irradiated by two laser beams. Molecular 
losses, induced by one of the two lasers, are suppressed 
when the second laser is tuned into resonance with a 
bound state (see Fig.[T|). Fitting a coupled-channel model 
to our experimental data, we have constructed an accu- 
rate Born-Oppenheimer potential for a 3 E+. This en- 
ables us to calculate the wave functions of triplet bound 
states as well as their binding energies to 60 MHzxh accu- 
racy (where h is Planck's constant) over the whole man- 
ifold of vibrational and low rotational (N<5) levels. Our 
theory and data agree to the extent that further refine- 
ment of the theoretical model requires a reduction in the 
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FIG. 1: Dark-state spectroscopy scheme for the 87 Rb2 a 3 E+ 
potential. Lasers 1 and 2 couple the molecular levels |i) and 
\v) to the excited level |e) with Rabi frequencies f2i,2, respec- 
tively. Laser 1 is kept on resonance while laser 2 can be tuned 
to any level of the o 3 E+ potential. Inset: Bound state level 
of Feshbach molecules as a function of magnetic field B z . The 
dashed line gives the position of the Feshbach resonance. The 
dot marks the Feshbach molecule state used in the experi- 
ments. 



experimental uncertainty. 

Molecular spectroscopy with cold atomic gases goes 
back to the beginnings of laser cooling |L2Hl4| . An ex- 
periment closely related to ours is the one by Araujo et 
at [1 51 ] where the Na 2 triplet ground state was explored 
using a magneto optical trap combined with two-color 
photoassociation spectroscopy. Our spectroscopy also 
makes use of a two-photon transition, but starts from 
cold Feshbach molecules rather than from free atoms. 
Recent investigations of the a 3 S+ potential of Rb 2 in- 
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elude the work of [16| . jrTj and [18|. Using one-color 
photoassociation of the a 3 S+ potential of laser cooled 
85 Rb, reference 16] put tight constraints on the position 
of the repulsive wall of the a 3 S^ potential but was not 
able to resolve the vibrational structure. The two other 
groups determined several ro- vibrational levels using flu- 
orescence spectroscopy fl7| or pump-probe photoioniza- 
tion spectroscopy of Rbo formed on helium nanodroplets 
18]. Our work goes well beyond these measurements as 
we fully resolve hyperfine, rotational and Zeeman struc- 
ture for almost all vibrational states. Highly precise data 
of the asymptotic behavior of the coupled a 3 S+- X 1 E+ 
system is contained in the large set of observed Feshbach 
resonances fl^l - E^ and in the two-photon photoassocia- 
tion measurements of four weakly bound levels at zero 
magnetic field by 23]. 

This article is organized as follows: Section [TTI presents 
the experimental setup and typical dark state spec- 
troscopy scans. Section IIIII discusses the relevant quan- 
tum numbers of our studies and the assignment of the 
observed lines. Section IIVI is a short summary of the 
coupled-channel model and the optimization procedure 
of the Born-Oppcnheimer potentials. SectionlVl discusses 
the progression of the substructure of the vibrational 
manifolds. We conclude the paper with a summary and 
an outlook towards further experiments in section IVTl 



II. EXPERIMENTAL SETUP AND DARK 
STATE SPECTROSCOPY 



The starting point for our experiments is a 50 /im- 
size pure ensemble of 3 x 10 4 weakly bound Feshbach 
molecules which have been produced from an atomic 
Bose-Einstein condensate of 87 Rb by ramping over a Fes- 
hbach resonance at a magnetic field of 1007.4 G (1G = 
10~ 4 T) [H[. They are trapped in the lowest Bloch band 
of a cubic 3D optical lattice with no more than a single 
molecule per lattice site (25[. The lattice depth for the 
Feshbach molecules is 60 £V, where E r = ir 2 fi 2 /2ma 2 is 
the recoil energy, with m the mass of the molecules and 
a = 415.22 nm the lattice period. Such deep lattices pre- 
vent the molecules from colliding with each other, which 
suppresses collisional decay. We observe lifetimes of a few 
hundred ms. After producing the Feshbach molecules, 
the magnetic field is set to 1005.8 G where the spec- 
troscopy is carried out. At this magnetic field, the bind- 
ing energy of the Feshbach molecules is 4.4(3) MHz x h 
(Fig. [T] inset). 

Our dark state spectroscopy works as follows: The Fes- 
hbach molecules in state |i) are irradiated by simultane- 
ous rectangular pulses from lasers 1 and 2. The pulses 
typically last 10 /is with Rabi frequencies f^i and Q2, re- 
spectively (Fig.[TJ. We keep laser 1 resonant with the 
|i) — |c) transition and at a power I\ of about 0.1 mW 
= 2ttx 0.3 MHz) such that in the absence of laser 
2 about half of the molecules are lost by spontaneous 
emission from |e). Laser 2 with its power I2 (up to a 
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FIG. 2: (a) Binding energies Eb(v) for the state a , where 
v is the vibrational quantum number. The line is the result 
of a coupled-channel model calculation after optimization of 
the Born-Oppenheimer potential. Five levels were not mea- 
sured, (b) Residues i. e. difference between experimental data 
and theory. Large error bars belong to early measurements 
without simultaneous wavemeter readings of both lasers. 



few hundred mW) is scanned. As long as laser 2 does 
not hit an |e) — \v) resonance, laser 1 will continue to 
induce losses. However, when an |e) — \v) resonance oc- 
curs, the initial state |i) is projected onto a dark state 

\^dark) = (0 2 |i) - ni|u))/v / Of + nf. Molecules in this 
dark state are shielded from excitation to the short-lived 
level |e) [2(|. This leads to a suppression of molecular 
losses. After the lasers are switched off we measure the 
number of molecules via a reverse magnetic field sweep 
through the Feshbach resonance, dissociating the remain- 
ing molecules in |i) into atoms which are detected by ab- 
sorption imaging. These measurements are destructive, 
and for each point in a scan, a fresh sample of Feshbach 
molecules has to be prepared. 

The level |e) has an excitation energy of about 295 
THzx/i with respect to |i) and a width T = 2ttx 8 MHz. 
Laser 1, a grating-stabilized diode laser, is Pound-Drever- 
Hall locked to a cavity which is in turn locked to an 
atomic 87 Rb-line. Laser 2, a TkSapphire laser, is free- 
running and typically drifts over a frequency range of a 
few MHz within seconds. Both lasers have a short-term 
laser line width of several tens of kHz. Their beams have 
a 1/e 2 intensity waist radius of 130 /im at the molecular 
sample, through which they propagate nearly collinearly. 
They are polarized parallel to the magnetic bias field B z 
(pointing in the vertical direction) and thus can only in- 
duce 7T transitions. The frequencies of both lasers are 
automatically read out using a commercial wavemeter 
(WS7 from HighFinesse) . The binding energy Eb of \v) 
minus the binding energy of our initial Feshbach state 
(4.4MHzx/i) corresponds to the frequency difference of 
the two lasers. The wavemeter has a nominal accuracy 
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of 60 MHz after calibration. Over the course of days we 
have observed drifts of ±200 MHz, e.g. by repeatedly 
addressing the same spectroscopic line. Over the length 
of a few experimental cycles (5 minutes) the wavemctcr 
is stable to within 10 MHz, which represents a random 
noise floor. Assuming a sufficiently smooth behavior of 
the wavemeter, drifts of the wavemeter typically affect 
the frequency measurements of laser 1 and laser 2 in a 
similar way, especially since the laser frequencies only 
differ by 3%. Thus, such common mode drifts cancel out 
in the binding energy to first approximation. Indeed, 
based on our experience with the wavemeter where we 
have measured binding energies of a few sharp lines over 
an extended period of time and via various intermediate 
levels we estimate the accuracy to reach 30 MHz. This 
includes peak position uncertainties due to variations of 
the number of molecules produced, as well as a frequency 
drift of laser 2 during the time between the laser pulse 
and wavelength measurement. 

Fig. [2] a) shows the measured binding energies of the 
triplet potential as a function of the vibrational quan- 
tum number v at a magnetic field of 1005.8 G. There are 
41 vibrational states with binding energies ranging from 
5 MHz x h to about 7038 GHz xh. The vibrational split- 
ting between the two lowest vibrational levels, v = and 
v = 1, is about 393 GHz. Each vibrational state has hy- 
perfine, rotational, and Zeeman substructure. This struc- 
ture is spread out over a range of about 20 GHz as shown 
in Fig. [3] for the states v = and 6. These spectra typi- 
cally consist of roughly 1000 points corresponding to an 
average step size of 20 MHz. Each point represents one 
production and measurement cycle which takes 28 s. 

In each spectrum of Fig. [3] we observe some 10 lines 
which vary markedly in line width. The width of each 
line is determined by the coupling between the levels |e) 
and \v), i.e. the Rabi frequency 0,2- Interestingly, for 
our measu ring scheme this width scales as f2| [23] and 
not as O2 [28l. [29I]. as one might expect. We have taken 
advantage of this enhanced broadening when searching 
for lines and vibrational manifolds, which otherwise can 
be like looking for a needle in a haystack. For example, 
for v — and an intensity I2 of a few hundred mW we 
reached line widths of several GHz. The substructure of 
the desired vibrational level then appears essentially as 
a single broad line with a width of about 20 GHz. Once 
this level was found the power was reduced, in order to 
resolve its substructure. 

In general, we expect spectra of different vibrational 
levels to be very similar. Up to v = 35 this is indeed 
the case. For v > 35 the vibrational manifolds start to 
overlap, as the splitting between them becomes smaller 
than 20 GHz. Spectra a) (v = 0) and b) (v = 6) of Fig.[3] 
are clearly similar. Some visible differences are artifacts, 
e.g. it seems that in spectrum b) two narrow lines are 
missing at Es/h s» 7021GHz. This can be explained by 
the fact that the lines were narrower (width < 10 MHz) 
than the local step size of that scan and thus escaped 
observation. 
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FIG. 3: Scans of two vibrational levels: a) v — 6, b) v = 0. 
Plotted is the remaining molecular fraction as a function of 
the binding energy Eb, which basically corresponds to the 
laser difference frequency. The scans were recorded using an 
excited state |e) with l g character. The labels s and d indicate 
rotational states N = and 2, respectively. The numbers 
after the labels s or d indicate the position in the spectrum. 



III. QUANTUM NUMBERS AND 
ASSIGNMENT 

A. Quantum numbers 

For a deeper understanding of the structure of the 
spectrum and its assignment we now discuss the relevant 
quantum numbers and selection rules. 

For the weakly bound levels like the Feshbach state 
or the vibrational levels close to the atomic asymptote 
Hund's coupling case (e) with atomic quantum numbers 
is most appropriate. Here none of the angular momenta 
couple to the molecular axis. The state vector is de- 
scribed by 

\(f A ,fn)f,l,F,M), 

where /a,/b are the total angular momenta for atoms 
A and B, / is the sum of both atomic angular momenta, 
I is the mechanical rotation of the atomic pair, F is the 
total angular momentum of the pair and M is its projec- 
tion onto a space fixed axis. At low magnetic fields the 
prepared Feshbach state |i) can be approximated by the 
state vector \ f A = 2, f B = 2, f = 2, 1 = 0, F = 2, M = 2). 
At the magnetic field used in the experiment (1005.8 G), 
F is no longer a good quantum number because of only 
weak coupling between / and /. A more appropriate state 
vector is 

|(/a,/b)/, m f ,l,mi,M). 

M is a good quantum number, while for example, /a and 
Jb have expectation values (3fJ of about 1.79 instead of 
quantum numbers /a = /b = 2 for G. 

Due to the strong hyperfine coupling of Rb and the 
large exchange energy, deeply bound levels of the triplet 
state can be described by Hund's coupling case (b^) at 
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low magnetic fields [3l| . namely: 

\N, (I,S)f,F,M), 

where / and S are the total nuclear and electronic spin 
quantum numbers, N the molecular rotation including 
electron orbital angular momentum, and F and M have 
exactly the same meaning as in Hund's case (e). The in- 
dex (3 indicates that the nuclear spin I is not coupled to 
the molecular axis. Since both atoms are in the 5S con- 
figuration of Rb, the molecular electronic orbital angular 
momentum is zero. This means that N — I and that / is 
the same as in Hund's case (e). Here we get 

\N,m N ,(I,S)f,m f ,M), 

as an appropriate state vector for higher magnetic fields, 
where M = tojv + to/. One can show that owing to 
the antisymmetry of the molecular wave function with 
respect to nuclear exchange (nuclear spin of 87 Rb i = |), 
molecules with even (odd) N in the a 3 E+ state must 
have either a total nuclear spin I = 1 or 3 (I = or 2) 
[il], HH- For the X 1 S+ singlet ground state this relation 
is reversed because it has g symmetry in contrast to the u 
symmetry of the triplet state. For large magnetic fields, 
/ loses its meaning as S and / start to decouple. 

Expectation values for the total nuclear spin and the 
electron spin for our Feshbach level are I — 1.56 and S = 
0.76, respectively. Thus we have significant electronic 
singlet-triplet mixing and consequently also mixing of the 
basis vectors with different I. In contrast, the excited 
intermediate state |e) has well defined quantum numbers 
S and /. Thus |e) largely determines which quantum 
numbers the deeply bound a 3 £^ levels will have in the 
Raman transition. 

The intermediate level |e) is located in the v' = 13 
manifold of the l 3 S+(5S , 1 / 2 + 5Pi/ 2 ) potential. Due to 
significant effective spin-spin interaction, the vibrational 
manifold is split into two components, l g and g , sep- 
arated by 47 GHz. As intermediate level |e) we cither 
choose a level with l g character or with 0~ character. 
Its rotational energy must be low because the Feshbach 
level has the lowest rotational quantum number N = 0. 
As rotation is low in |e), decoupling of S and / from the 
molecular axis by rotation is not yet important. We can 
then approximate the state |e) by quantum numbers for 
Hund's case (a Q ) where the index a indicates that the 
nuclear spin I is coupled to the molecular axis, 

\SE,ISli,F,M). 

The projections of the electronic and nuclear spin onto 
the molecular axis appear as quantum numbers E and 
fli, respectively. 

The |e) level with l g character is energetically the low- 
est within the l g (v = 13) manifold. It has the quantum 
numbers S = 1, |£| = 1, / = 3, « 3, F 2 and 
M = 2 for low magnetic fields B z . Its excitation energy 
is 294.6264(2) THzx/i with respect to |i) at a field of 
1005.8G. 



The |e) level with g character has an excitation en- 
ergy of 294.6736(2) THz. It can be described by the quan- 
tum numbers S — 1, £ = 0, I = 3, M = 2 and F = 3 at 
B z = 0. Because of its low hyperfine coupling the better 
choice of basis vector here is 

|5E, (JI)F,M), 

where J results from the coupling of the electronic spin 
and the molecular rotation (J = N + S). For our 0~ 
|e) level J is approximately zero. Our l g |e) state is a 
superposition of N=l and 3. 

As stated before, due to the laser polarization along 
the magnetic field only tt transitions are allowed, which 
results in the selection rule AM = 0. Further, in a one- 
photon transition parity has to change. The Feshbach 
state |i) has a total parity "plus" because it is a E + state 
and (— 1) N = 1. Thus we can only address |e) levels with 
"minus" parity and \v) levels with "plus" parity. This 
means that the \v) level must have an even rotational 
quantum number TV = 0, 2, 4, .... For the |e) level only 
quantum numbers N = 1, 3, ... or superpositions of these 
are available. In fact, as the 1 3 E+ state is well described 
in a Hund's case (a) basis, N is in general not a good 
quantum number for the |e) level. The selection rule 
AiV = ±1 for N determines the range of reachable levels 
for \v) according to the superposition in |e). The selection 
rules AI = and AS — are important for the transition 
|e) to \v). For the transition |i) to |e) they are, however, 
nearly irrelevant since I and S are not good quantum 
numbers for level |i). 

B. Assignment of spectral lines 

Figured] shows measured and calculated (based on the 
coupled-channel model, see section HV]) lines of the v = 6 
spectrum, where the excited state |e) with 0~ character 
was selected. The lines form three groups with the quan- 
tum numbers / = 4, 3, and 2, according to the hyperfine 
coupling of S = 1 and 1 = 3. This is a clear indication 
that the hyperfine energy is still dominant compared to 
the Zccman energy. 

Each group consists of one line with N = corre- 
sponding to a non-rotating molecule, and several lines 
with N — 2 which are shifted to lower binding energy 
by about 2 GHz x h due to rotation. The fact that we 
do not observe lines with N > 2 can be explained as 
follows: The excited state |e) has "minus" total parity. 
Since J = and S = 1, N must be equal to one and is 
also a good quantum number for this lowest level in the 
0~ manifold. Thus when using the state with g charac- 
ter only final states \v) with N = 0, 2 can be addressed 
in the a 3 E+ state due to the selection rule AN = ±1. 

The overall structure of each vibrational level can be 
understood with a relatively simple effective Hamiltonian 

H=^S-I + B V N 2 + fig g s S z B z 
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FIG. 4: Scan of the hyperfine, rotational and Zeeman struc- 
ture in the a 3 E^(w = 6) manifold using an excited level |e) 
with OJ -character. The x-axis shows the binding energy Eb- 
The quantum number / is shown below each of the three 
groups of lines. Thick black lines above the spectrum indi- 
cate states with N = 0, grey lines correspond to states with 
N = 2. The upper row of numbers indicates the expectation 
values of the magnetic quantum number m/ at the magnetic 
field of 1005.8 G. 



with the atomic hyperfine structure constant A — 3.42 
GHzx47r 2 //i [34J, the total electronic and nuclear spin op- 
erator S and I, respectively, and the operator for molec- 
ular rotation N. B v is the rotational constant of the 
desired vibrational level v. The last term describes the 
Zeeman effect of the electronic spin when exposed to an 
external magnetic field B z in z-direction. The nuclear 
Zeeman term is neglected. The Zeeman effect can be 
evaluated for the case of strong hyperfine coupling such 
that / (with f = S + I) is still a good quantum number. 
Using a simple vector model for f one can derive a Lande 
factor gf. 
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For each /, the N = 2 group is shifted to lower bind- 
ing energy by 6B V compared to N = 0, corresponding to 
about 2 GHz for the v = 6 manifold (see Figure 01 . Each 
N = 2 group is split by the Zeeman energy according to 
the m,f quantum number of each line. Only those lines 
which have quantum numbers where M = 2 = rn t + mjv 
can be observed. The level / = 3 has a small gf factor 
(1/6 for 1 = 3), which gives rise to a small Zeeman split- 
ting. Additionally, one sees that the splitting between 
/ = 4 and / = 2, which would be about 12 GHz for the 
pure hyperfine part, is enlarged by a Zeeman contribu- 
tion because of the different signs of gf for / = 4 and 
f = 2. In Figure [4] the expectation values for rrif re- 
sulting from the coupled channel model are given in the 
upper row and show that these lead to good quantum 
numbers for / = 4 and / = 2 but are less good for / = 3, 



where the hyperfine and electronic Zeeman energies are 
small. The mixing of states with different nif for / = 3 
is not included in the simple model of Eq. Q] 

As discussed before, for the intermediate state with 0~ 
character we exclusively observed N = and N = 2 lev- 
els in the a 3 E+ potential. This restriction does not nec- 
essarily apply when using the intermediate level |e) with 
l g character. The corresponding spectrum in Figure [3] a) 
shows the same v — 6 manifold as in Figure 21 only |e) 
has l g instead of 0~ character. It has additional lines 
(at ~4838 GHz and at ~4832 GHz) which match the pre- 
dicted positions of N = 4 levels. 

As mentioned in Section llll Al the level |e) with l g char- 
acter would be described by a superposition of N = 1 and 
N = 3 states in a Hund's case (b) basis. We thus expect 
to see transitions to levels with A = 0, 2, 4 of the a 3 E+ 
state according to the selection rule AN = ±1. In order 
to avoid confusion, we note that the N = and N = 2 
lines between 4825 GHz and 4830 GHz which are clearly 
visible in Figure 2] are weak or are not even observed in 
Figure |3] For the same reasons as discussed at the end 
of section |lTl the experimental step size of about 10 MHz 
might have been too large compared to the narrow line 
widths for these transitions to be seen. 

It was not always necessary to carry out a complete 
scan as in Figure 0] in order to assign quantum numbers 
to observed lines in arbitrary vibrational levels. Often it 
was sufficient to measure a few characteristic lines and 
splittings and to compare them to the calculated spec- 
trum. These data were then used to optimize the coupled 
channel model along with its a 3 S+ and A 1 !]^ Born- 
Oppenheimer potentials. 



IV. COUPLED-CHANNEL MODEL AND 
OPTIMIZATION OF THE a 3 E+ POTENTIAL 



In our work we use a coupled-channel model |35M3 
that can calculate all bound states of the X 1 £+ and 
a 3 E+ states, which correlate with the atomic asymptote 

(1) 

5 2 5i/2 + 5 2 5i/2- The program has helped us with the 
search for lines as well as with their identification. Using 
our data we are able to optimize the Born-Oppenheimer 
potential of the a 3 E+ state as well as to improve the 
potential of the ground state X given in [30|. In the 
following we briefly describe the model and explain how 
the A 1 S+ and a 3 S^ Born-Oppenheimer potentials are 
optimized. 

To cover the full range of experimental data by a single 
theoretical model, only a coupled channel analysis is ade- 
quate. It includes the calculation of the molecular bound 
states as well as the scattering resonances. It takes into 
account the A 1 !]^" and a 3 S+ potential functions, the 
hyperfine coupling, the Zeeman interaction, rotation and 
the effective spin-spin interaction. Such a theoretical ap- 
proach is described in several papers, e.g. (38[. 

For the present analysis we include our measurements 
covering 135 lines of the a 3 E+ state as well as data from 
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other work. For the X 1 £+ ground state we added data 
from Fourier transform spectroscopy by Seto et al. [39| 
with more than 12000 lines. We also include measure- 
ments of Feshbach resonances for the three isotopologues 
85 Rb 2 , 87 Rb 2 , and 85 Rb 87 Rb [H-El], the four asymptotic 
levels from |23j |. and measurements from Fourier trans- 
form spectroscopy for the a 3 £+ state reported by Beser 
et al. M- 

The full Hamiltonian (cf. |H, EH 53) for a pair of 
atoms A and B can be written in the form 

H = T n + U x (R)Px + U a {R){l - P x ) 
+ a A (R)sA ■ iA + a B {R)s B • %b 

+ [{gsASzA - giAlzA) + (g.sBS z B - 9iBizB)] ^bB z 



(2) 



In the present case of the homonuclear molecule 87 Rb2 
the parameters with index A are equal to those of index 
B. The first term in the first line shows the kinetic energy 
T n where we take the atomic masses from recent tables 
by G. Audi et al. [42| . The next two terms describe the 
potential energies Ux and U a for the motion of the atoms, 
where Px and 1 — Px are projection operators on to the 
uncoupled states X and a, respectively. R corresponds to 
the internuclear separation of the two atoms. The second 
line shows the hyperfine interaction between the atomic 
electron spins sa,b and the atomic nuclear spins ia.b- 
The main contribution to the functions aA,B{R) is the 
Fermi contact term. The R dependence of the hyperfine 
parameters accounts for several effects: It takes into ac- 
count the electronic distortions of one atom by the other, 
i.e. the binding, and an effective coupling of the electron 
spin of one atom with the nuclear spin of the other atom. 
We start with i?-independent atomic couplin g co nstants 
taken from a compilation by Arimondo et al. [43| ■ These 
constants are later refined by a simple Ansatz for the R 
dependence (discussed later in Eq. [8]). We neglect the 
nuclear quadrupole moment in the hyperfine interaction, 
which might come into play for deeply bound levels. The 
third line in Eq. [2] gives the Zeeman energy from the 
coupling of the electron spin and the nuclear spin with 
an external homogeneous magnetic field B z in z direc- 
tion. The electronic and nuclear g-factors for the atomic 
ground state of the Rb isotopes are taken from the re- 
port in [4§ |. This term couples states with different / 
quantum numbers. The last line contains the spin-spin 
interaction which couples different N states of basis (b). 
It is formed by the total molecular spin S and its pro- 
jection on the molecule fixed axis Z . The parameter A is 
a function of R, of which one part has 1/i? 3 dependence 
as a result from the magnetic dipole-dipole interaction. 
In addition, A contains contributions from second order 
spin-orbit interactions. The final analysis showed that 
such a contribution is significant within the achieved ex- 
perimental accuracy. For example, this was important 
for the precise location of Feshbach resonances involving 
/ = 1 and / = 2 bound levels. 



The functional form of the two Born-Oppenheimer po- 
tentials, X and a 3 £+, is split into three regions on 
the internuclear separation axis R: the short range repul- 
sive wall (R < Rsr), the asymptotic long range region 
(R > Rlr), and the intermediate deeply bound region 
in between. The analytic form of the potentials in the 
intermediate range, Uir, is described by a finite power 
expansion of a nonlinear function £ which depends on the 
internuclear separation R, 



i=0 



R — Rm 

R + bRr, 



(3) 



(4) 



Here the Oj are fitting parameters (see table |TT|. We 
choose b and R m such that only few parameters Oj are 
needed for describing the steep slope at the short inter- 
nuclear separation side and the much smaller slope at 
the large R side. R m is chosen close to the value of the 
equilibrium separation. The potential is extrapolated for 
R < Rsr by the short range part Usr with: 



U S r(R) = mi + u 2 /R 



(5) 



We adjust the parameters u\ and u 2 to get a continuous 
transition at -Rsr. The final fit uses N s s» 4.5 for both 
the X 1 ^ and a 3 £+ states. 

For large internuclear distances (R > Rlr) we adopted 
the standard long range form of molecular potentials 



U L r(R) = -Ce/R 6 - Cs/R 8 - C w /R 10 ± E cxch , 
where the exchange contribution given by [44j 
E cxch = A ex R7 exp(~/3R) 



(6) 



(7) 



is negative for the singlet and positive for the triplet po- 
tential. By adjusting the parameter ao in Eq. [3] we can 
assure a continuous transition from C/lr to £/ir. As men- 
tioned, the data on hand include three different isotopo- 
logues, namely 85 Rb 2 , 87 Rb 2 , and 85 Rb 87 Rb. Using a 



model developed earlier for LiK 45] we checked in the 
final calculations that the data are not sufficiently pre- 
cise to extract deviations from the Born-Oppenheimer 
approximation. Thus we derive the potentials without 
any mass scaling correction. 

Having discussed the complete physical model, we are 
now ready to calculate all the relevant bound state en- 
ergies and scattering properties to compare them with 
the experimental data. We have decided to evaluate the 
Hamiltonian in Hund's basis (e) as / is still a relatively 
good quantum number. This is due to the fact, that the 
hyperfine interaction is larger than the Zeeman interac- 
tion in our experiment. As the total electron spin is not 
a good quantum number in Hund's case (e), the choice 
of this basis leads to significant non-diagonal matrix ele- 
ments from the Born-Oppenheimer potentials, which are 
given for a pure singlet or triplet state. 



7 



We evaluate the free parameters of the model with a 
self-consistent iteration loop and alternate between: (i) 
coupled channel calculations of the full Hamiltonian and 
(ii) solving the Schrodinger equation separately for the 
states X 1 S+ and a 3 S+ using only the first line of Eq.[2] 
and applying the Numerov procedure. The coupled chan- 
nel calculations in step (i) are used to determine the hy- 
perhne and Zeeman structure. From this we construct 
hyperfme free spectroscopic data. This data is the in- 
put for the Born-Oppenheimer potential fits in step (ii). 
We incorporate the fitting routine for the plain Born- 
Oppenheimer potentials in step (ii) . The optimization of 
the singlet and triplet potentials is done simultaneously 
as both potentials have a common asymptote. These 
asymptotic potentials are given in Eq.|6] with equal dis- 
persion coefficients. 

Normally, the iteration loop between the potential 
function fit and the coupled channel calculation for pro- 
ducing hyperfme structure free data converges in a few 
steps, but we observed some systematic deviations be- 
tween measurements and calculations. Specifically the 
hyperfme splitting showed variations of a few percent 
within the vibrational ladder. Even though our exper- 
imental accuracy is not better than 30 MHz such small 
variations are observable since the hyperfme splitting in 
87 Rb 2 of 12 GHz is so large. 

In order to theoretically account for the hyperfme vari- 
ations we extended the model. We changed the fixed 
atomic hyperfme parameters to a function in R. We 
chose a function which switches at a distance Rq from 
the atomic value of the hyperfme constant to another 
value for a deeply bound dimer, 

a A /B = a Rb ^1 + e(R _ Ro C } } /AR + 1 ^ ■ (8) 

Here a Rb is the atomic hyperfme constant, c/ the frac- 
tional change of the constant and Ro and AR describe the 
switching distance and its width, respectively. We also 
tried several other simple switching functions, which pro- 
duced about equal fit quality. The function in Eq.[S]is eas- 
ily applicable for other isotopes by introducing the proper 
atomic hyperfme constant, because the scaling c/ will be 
independent of the isotope. We chose Ro = 11.0 ao and 
AR, = 0.5 a (where a = 0.5292 x 10" 10 m is the Bohr 
radius), such that switching takes place approximately 
at the minimum of the a 3 X+ potential. From our fits we 
obtained an amplitude c/ = —0.0778, which corresponds 
to a variation of the hyperfme coupling across the poten- 
tial depth of up to 8%. The influence of this hyperfme 
variation within an individual vibrational manifold will 
be smaller because of the averaging of Eq.[8] over the vi- 
brational wave function. As mentioned before, also the 
spin-spin interaction needs optimization in order to ex- 
plain systematic shifts of Feshbach resonances in s-wave 
and p-wave scattering channels, as measured in [l9j . The 
spin-spin interaction couples different partial waves I sub- 
ject to the selection rule Al — 0, ±2 such that resonances 
in a s-wave scattering channel involve bound states with 



I = and 1 = 2. The spin-spin interaction splits the res- 
onances according to |m;|, the projection quantum num- 
ber of rotation on the space fixed axis. We use a simple 
functional form for A(i?) in Eq. [5] 

A(i?) - - '^a 2 f-L + a SO exp {-b{R - R so ))\ , (9) 

that consists of two terms. The first term represents the 
magnetic dipole-dipole interaction and the second term is 
the second order spin-orbit contribution. If Eq.[9]is given 
in atomic units, a is the fine structure constant. Since 
the few data at hand cannot be highly sensitive to the 
actual function, we adopted values for b and Rso from 
a theoretical approach by Mies et al. for Rb 2 |46| (b= 
0.7196 a^ 1 and Rso= 7.5 ao). We fitted the parameter 
aso, which yielded aso— -0.0416 &q 3 . With these param- 
eters the second part in Eq. [9] contributes significantly to 
the effective spin-spin interaction in the internuclear sep- 
aration interval R < 20ao. This affects bound levels of 
the triplet state that determine the Feshbach resonances. 
After optimization we achieve an accuracy of 0.1 G for 
the position of the Feshbach resonance. This correction 
to the constant A does not influence the description of 
the other bound states within their uncertainty. 

Finally, we found that we can improve the fit by adding 
to the long range formula Eq. 6 a term of the form 
—C26/R 26 , and retaining all constraints such as conti- 
nuity as before. From our fits we obtain an amplitude 
C26 which contributes about a thousandth of the total 
long range energy at the connection point Rlr- The ex- 
ponent of 26 was chosen so that the term is negligible 
outside a small region around the long range connection 
point R = Rlr- 

The final parameter sets of the potentials are shown 
in table U and |H] for the singlet and triplet states, re- 
spectively. The derived potentials and the corrections 
defined above agree very closely with all observations to 
within their experimental uncertainties, and the normal- 
ized standard deviation (i.e. standard deviation / exper- 
imental uncertainty) is close to one. In fact, depending 
on which of our experimental data sets (hyperfme free 
spectra, binding energies from dark state spectroscopy 
or Feshbach resonances) we compare to the model cal- 
culations, the normalized standard deviations only vary 
slightly, ranging from 1.01 to 1.3. The normalized stan- 
dard deviation from the joint calculation over the huge 
body of ro-vibrational energies (12459 data points from 
the states X and a 3 S^) is about 1.01, which is quite 
satisfactory. 

As another result of our analysis we were able to elim- 
inate an ambiguity in the rotational assignment of the 
Fourier transform spectroscopy data reported by Beser 
et al. for the state a 3 S+ [17]. The authors stated in 
the rotational assignment an ambiguity by A A = ±1. 
We determined that the shift of the rotational quantum 
number must be AN = +1. Afterwards we used the data 
from (l7| with that assignment for the further fits. The 
result in Table [TTI includes these data and will give a bet- 
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ter confidence in predicting rotational levels with higher 
rotational quantum numbers. 

V. PROGRESSION OF VIBRATIONAL LEVELS 
AND THEIR SUBSTRUCTURE 

In the following we will discuss interesting insights 
which we have gained from our analysis of the coupled 
system. In particular we investigate the progression of 
several quantities in the vibrational ladder and the mix- 
ing of singlet and triplet states. We concentrate on the 
details of the a 3 S+ state, because it was studied in full 
resolution of hyperfine and Zeeman energy over the whole 
vibrational ladder for the first time. Such a body of data 
does not exist for any other alkali dimer. 

A. Vibrational ladder and rotational progression 

We return to the vibrational states shown in Fig. [2] 
which allows us to compare in detail the optimized 
coupled-channel model with our experimental findings for 
all vibrational states. 




U ' 1 1 1 1 1 1 

5 10 15 20 25 30 35 

Vibrational quantum number v 

FIG. 5: Energy splitting AE B between the "sl"-level (JV = 0) 
and the "dl"-level (N = 2) for a given vibrational quantum 
number v. Large error bars correspond to early measure- 
ments without simultaneous wavemeter readings (see text). 
The continuous line is a calculation based on the coupled 
channel model. 

The binding energies given in Fig. [2] a) correspond to 
the most deeply bound level in each vibrational level, i.e. 
the state with quantum numbers N = 0, / = 2 (at G) 
and M = 2. We will refer to this level as "si" . The "si" 
level of v = in the a 3 £+ potential is also the lowest 
bound state in that potential and has an observed bind- 
ing energy of (7038.067±0.050)GHzx/i at 1005.8G with 
respect to the lowest atomic asymptote wia — 1 

and fs ~ 1 , itlb = 1 • A list of all the other measured 
and calculated bound state energies together with their 
expectation values of the quantum numbers can be found 
in the supplementary material (47| . 

Fig.[5]b) shows the residues between our experimental 
"si" data and the optimized model. In general the model 



TABLE I: Parameters of the analytic representation of the 
potential of state X £j~. The energy reference is the disso- 
ciation asymptote and the term energy is the depth of 
the potential. Parameters with * are set for continuous ex- 
trapolation of the potential and those with ** see text for 
interpretation. 



R < R S r = 3.126 A 


N s 




4.53389 


u l 




-0. 638904880x10* cm" 1 






0.112005361X10 7 cm -1 A JV « 




Rsk 


< R < i?, L R = 11.000 A 


b 




-0.13 


Rm 




4.209912760 A 


ao 




-3993.592873 cm" 1 


«i 




0.000000000000000000 cm" 1 


a ■> 




0.282069372972346137xl0 5 cm" 1 


03 




0.560425000209256905 xlO 4 cm" 1 


a 4 




-0.423962138510562945 x10 s cm" 1 


05 




-0.598558066508841584x10 s cm" 1 


06 




-0.162613532034769596 x10 s cm" 1 


07 




-0.405142102246254944x10 s cm" 1 


08 




0.195237415352729586x10 s cm -1 


ay 




0.413823663033582852x10 s cm" 1 


O10 




-0.425543284828921501 xlO 7 cm" 1 


an 




0.546674790157210198x10 s cm" 1 


012 




0.663194778861331940x10 s cm" 1 


ai3 




-0.558341849704095051x10 s cm" 1 


014 




-0.573987344918535471 xlO 9 cm" 1 


ai5 




0.102010964189156187X10 10 cm" 1 


«16 




0.300040150506311035 xlO 10 cm" 1 


«17 




-0.893187252759830856 xlO 10 cm" 1 


«18 




-0.736002541483347511 xlO 10 cm" 1 


ai9 




0.423130460980355225 xlO 11 cm" 1 


«20 




-0.786351477693491840X10 10 cm" 1 


«21 




-0.102470557344862152 xlO 12 cm" 1 


022 




0.895155811349267578X10 11 cm" 1 


«23 




0.830355322355692902 xlO 11 cm" 1 


a24 




-0.150102297761234375 xlO 12 cm" 1 


<J25 




0.586778574293387070X10 11 cm" 1 


R > -Rlr 


c 6 




0.2270032x10 s cm" 1 ! 6 


c 8 




0.7782886x10 s cm^A 8 


Cio 




0.2868869X10 11 cm" 1 ! 10 


^26 




0.2819810xl0 26 cm-iA 26 


A ex 




0.1317786x10 s cm^A-T 


7 




5.317689 


P 




2.093816 A" 1 


Derived constants: 


equilibrium dist 


ance: 


R*= 4.20991(5) A 


electronic term 


energy 


T e x = -3993.5928(30) cm" 1 
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TABLE II: Parameters of the analytic representation of the 
potential of state a 3 E+. The energy reference is the disso- 
ciation asymptote and the term energy T" is the depth of 
the potential. Parameters with * are set for continuous ex- 
trapolation of the potential and those with ** see text for 
interpretation. 



R < R SR = 5.07 A 


N s 




4.5338950 


u l 




-0.619088543xl0 3 cm -1 


u* 




0.956231677xl0 6 cm~ 1 A Ns 




Rsn 


< R < i? L R = 11-00 A 


b 




-0.33 


Km 




6.0933451 A 






-241.503352 cm' 1 


"1 




-0.672503402304666542 cm" 1 


U-2 




0. 195494577140503543 xlO 4 cm' 1 


Q3 




-0. 141544168453406223 xlO 4 cm" 1 


a 4 




-0.221166468149940465X10 4 cm" 1 


a 5 




0.165443726445793004 xlO 4 cm" 1 


<Mi 




-0.596412188910614259 xlO 4 cm" 1 


"7 




0.654481694231538040 xlO 4 cm^ 1 


as 




0.261413416681972012x10 s cm -1 


a<) 




-0.349701859112702878x10 s cm" 1 


OlO 




-0.328185277155018630x10 s cm' 1 


an 




0.790208849885562522x10 s cm -1 






-0.398783520249289213x10 s cm -1 


R> H LR 


c 6 




0.2270032xl0 8 cm" 1 ! 6 


c 8 




0.7782886xl0 9 cm' 1 A 8 


Cio 




0.2868869X10 11 cm^A 10 


^26 




0.2819810X10 26 cm" 1 ! 26 






0.1317786x10 s cm _1 A~ 7 


7 




5.317689 


P 




2.093816 A -1 


Derived constants: 


equilibrium distance: 


R%= 6.0940(10) A 


electronic term 


energy 


T e a = -241.5034(30) cm" 1 



agrees very well with the measurements to within the er- 
ror bars. The data points with larger error bars belong 
to early measurements without simultaneous wavemeter 
reading of both lasers, which leads to a significant in- 
crease in the experimental uncertainty (see section Hi]) . 

We now investigate the progression of the rotational 
splitting in the vibrational ladder. For this, we consider 
the "si" level (N — 0) and its nearest neighbor in our 
spectra with (TV = 2) which we call "dl" (see Figs. [3] 
and Ul. The "si" and "dl" levels both have / = 2 and 
rrif = 2. Thus the splitting between them is to a high 
degree rotational energy. The splitting decreases with 
increasing v, which is due to the fact, that the mean dis- 
tance between the Rb2 nuclei (and hence the effective 
moment of inertia) increases as v increases. We have di- 
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FIG. 6: a) Progression of the substructure of the vibrational 
manifold at B — 1005.8 G. Shown are calculated a 3 Ej levels 
with quantum numbers J~3, M = 2, N = 0,2 and v = 
0... 34. Thick black lines correspond to N = 0, grey ones to 
N — 2. The "s2" level serves as energy reference (AEb ~ 0) 
in each vibrational substructure. In addition to the triplet 
levels nearby singlet levels of the X Ejj" potential are also 
shown (Thick black x's: N = 0, grey x's: N = 2). 
b) Symbols as in part a). To discuss the singlet-triplet mixing 
from a theoretical point of view we now also include lines 
with I « 1 but only iV = for clarity. On the bottom, the 
approximate quantum numbers (J, /) are given. See text for 
details. 



rectly observed this behavior in our experiments (see Fig. 
[5]). There is very good agreement between the experimen- 
tal data and the calculation (continuous line) using the 
optimized potential of the a 3 E+ state. 

Besides the "si" level (N = 0, / = 2) and the "dl" 
level (N = 2, / = 2) we also observe states with TV = 4 
and / — 2 for several low lying vibrational levels. The 
two lines in Fig. 3 a) at 4839 GHz and the two lines in 
Fig. 3 b) at 7032 GHz have N = 4. Observation of these 
levels improved the precision when fixing the position of 
the a 3 E J potential minimum in terms of the internuclear 
distance, or the effective rotational constant B v . This 
turned out to be important for the reassignment of the 
observations by Beser et al. [17], as discussed at the end 
of section IV. 
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B. Hyperfine splitting and singlet-triplet mixing 

It is instructive to investigate the progression of the 
multiplet structure of the vibrational levels. Figure [6]a) 
shows calculated levels for M = 2 stacked on top of each 
other for increasing v of the a 3 E+ potential. We concen- 
trate on the triplet states (S = 1) with rotation N = 0, 2 
and we restrict ourselves further to triplet levels with 
1 = 3. Additionally, we show singlet levels (S = 0) with 
1 = 2 that are located in the vicinity of the triplet lev- 
els. Thus, the typical stick spectrum of each vibrational 
manifold in the a 3 E+ state looks like the one in Fig. |4] 
There are three groups of lines with / = 4, 3, 2, respec- 
tively. The vibrational quantum number v runs from to 
34, where the multiplets do not yet overlap for different 
vibrational levels. 

In order to properly stack the data, we have chosen 
the level with N = 0, / ss 3 to be the energy reference 
(AE B = 0) for each vibrational level. We call this level 
"s2" (see also Figs. 3 and 4). This level is quite insensi- 
tive to mixing with singlet levels which makes it a good 
reference because there is no singlet level with / = 3 and 
even parity for direct coupling. It is clear from Fig. 6 
a), that the multiplet structure of different vibrational 
manifolds is similar, at least from v = to about v = 30 
as expected from the simple model Hamiltonian in Eq.[T] 
The structure changes in a smooth and monotonical way 
with v. For each of the quantum numbers / = 4, 3, 2 
the splitting between N = and N = 2 decreases with 
increasing v as discussed previously in section lV Al 

Appreciable mixing of a singlet and a triplet level can 
occur when the two levels with / = 2 are located energet- 
ically closely enough, in our case within a few GHz, cor- 
responding to the strength of the hyperfine interaction. 
One effect of mixing is a shift of the level positions due to 
level repulsion, clearly seen in Fig. [6] For v = 24, for ex- 
ample, the triplet lines of / = 2 are pushed to the left by 
about 0.4 GHz by the close singlet levels on the right hand 
side. For v < 20 we also observe narrow coincidences be- 
tween triplet and singlet levels, e.g. for v = 15. Here 
almost no perturbation appears in the graph, in contrast 
to the case v = 24. We attribute this to a significantly 
lower overlap of the vibrational wavefunctions of the sin- 
glet and triplet levels for the case v = 15 compared to 
that of v = 24. For vibrational quantum numbers v > 30 
mixing is very strong and happens for every vibrational 
level. Because the long range behaviour of the a 3 E+ and 
X states is similar, the overlap of the wave functions 
will become large for high v. The vibrational spacing will 
become similar to the hyperfine splitting, and the state 
vectors here are best described by Hund's coupling case 
(e) i.e. quantum numbers of atom pairs. 

Singlet-triplet mixing not only occurs for triplet 
molecules with 1 = 3 but also for I = 1. Figure H] b) 
shows (N = 0) triplet levels with 1 = 3 and 1 = 1 for 
v = 25 to 34. The repulsion of the 1 = 1 levels from 
the singlet levels is clearly visible. The figure shows an 
avoided-crossing-like behavior for the levels on the right 
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5 10 15 20 25 30 35 40 
Vibrational quantum number v 

FIG. 7: "Discrete second derivative" of Eb(v), i.e. 
S 2 E B (v) = E B (v + 1) - 2E B (v) + E B (v - 1). The solid curve 
connects the points from the "si" levels while the diamonds 
correspond to the "s2"-state. 

as a function of v, indicating the strong mixing between 
1 = 3 and 1 = 2 which also means that u/g symmetry 
is broken for these levels. Further, our calculations show 
that only levels with the same / and N quantum num- 
bers mix considerably. The singlet lines shown here have 
/ = 2 and N = 0. Apparently, despite the relatively 
strong magnetic fields of about 1000 G / is still quite a 
good quantum number. Indeed, / appears as quantum 
number in both state vectors for Hund's case (b) and 
(e) and loses its meaning only for much higher magnetic 
fields. N is good due to the small effective spin-spin in- 
teraction. 

The shift in position due to mixing can be traced more 
clearly with a simple difference method discussed in the 
following. Fig. [7] shows the "discrete second derivative" 
S 2 E B (v) of the function for the binding energies E B (v), 
i.e. S 2 E B (v) = E B {v + 1) - 2E B {v) + E B (v - 1). In 
other words it is the difference between the neighbor- 
ing energy splittings. The curve calculated from binding 
energies "si" exhibits sudden jumps for particular vibra- 
tional levels. These are due to singlet-triplet mixing and 
the v positions are consistent with Fig. [5] In contrast, 
the 6 2 E B (v) curve for the "s2"-state is smooth, and thus 
does not indicate mixing with the singlet lines, which 
justifies its choice as energy reference in Fig. [5] 

We have confirmed the singlet-triplet mixing experi- 
mentally. Fig. [5] shows scans of parts of the vibrational 
levels v = 28, 31, 33 where the intermediate level |e) with 
l g character was used (see also Fig. [3]). The "s2" -level is 
chosen to be the energy reference at AE B = as before. 
For v = 28 we observe a line structure similar to that 
of Figs. [3] and 2J Thus we take the v = 28 spectrum 
as reference of the pure case "triplet". In fact, from our 
calculations we see that the next singlet level is located 
about 30 GHz away. This detuning gives rise to only a 
very small mixing which slightly lowers the triplet char- 
acter of the "si" -level to 0.99. The singlet level obtains 
a triplet character of 0.02. However, due to the selection 
rule AS = 0, the transition to the intermediate state |e) 
with its very pure triplet character would simply be too 
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FIG. 8: Observation of strongly perturbed singlet levels. We 
compare sections of line spectra for the vibrational levels 
v = 28,31 and 33. As in Fig. [6] the "s2"-level is chosen as 
energy reference AEb = 0. The lines (and crosses) above the 
experimentally observed spectrum are from coupled-channel 
calculations. Black lines (crosses) represent N = levels, 
while thin grey lines (crosses) correspond to N = 2 levels for 
S = 1 (S = respectively). The lines that originate from the 
singlet states are labeled with " 1 S ". 



40 



JJ 


|J> 11'- 








d1 


CM 

v 30 


JJ 


s2 


I 1 

I II s 1 


CM 

E 
o 


<J2 




dl 


m 20 



- 10 

CM 



C. Franck-Condon overlap 















- 


i 3 sV 

g 


= 13)^a 3 2; + 

u 





























5 10 15 20 25 30 
Vibrational quantum number V 



35 



40 



FIG. 9: Normalized transition matrix element Q^/i^^/li) 
between the excited level |e) («' = 13, l g character) in the 
1 3 E+ potential and the "si" level with vibrational quantum 
number v in the lowest triplet potential a 3 E^. 



weak to see. 

This situation changes drastically for v = 31 and 33. 
Two additional lines are visible in the spectrum which 
originate from singlet states with M = 2, N = and 2. 
Here the singlet-triplet mixing is close to 40% which 
makes the singlet lines easily detectable. Additionally, 
for v = 31 the "si" and "dl" components are shifted to 
lower values than expected from the reference spectrum 
v = 28 due to repulsion by the singlet component on the 
high energy side. In the case of v — 33 the singlet compo- 
nent pushes the "si" and "dl" in the opposite direction. 
Our calculations show that even the chosen energy refer- 
ence "s2" starts to show mixing, indicated by its reduced 
triplet character of 0.98. In parallel to the singlet-triplet 
mixing, the quantum number I loses its meaning as well, 
e.g. the level "si" of v = 33 has an expectation value for 
/ of 2.75 instead of 3, due to a significant contribution of 
1 = 2 from the singlet state. Comparing the stick spec- 
trum line positions with observed lines in Fig. [5] we note 
the very nice agreement between experiment and theory. 
We want to mention that we measured the singlet levels 
at the predicted positions from first calculations. This 
emphasizes the predictive power of the presented model. 

The fact that we observe singlet-triplet mixing for rel- 
atively deeply bound levels with binding energies of a few 
hundred GHz x h is already quite interesting. In addition, 
it provides valuable information for fixing the energy po- 
sition of the triplet levels with respect to the singlets with 
high precision. This is especially important for the large 
body of data of the singlet system, see section IIV1 which 
was obtained with no connection to the triplet state [39j j . 



When scanning over all vibrational levels, the transi- 
tion matrix elements for the transition from |e) to \v) are 
not constant but oscillate as a function of v. This oscil- 
lation is mainly due to variations in the Franck-Condon 
overlap between the |e), v' = 13 vibrational wave func- 
tion and the vibrational wave functions \v) of the a 3 S+ 
potential. Figure [9] shows the normalized transition ma- 
trix element, C2 = ^li{v) / '(2ir yfh) , where 0,2 is the Rabi 
frequency and I2 is the intensity of laser 2. is deter- 
mined from the measured width of the dark resonance 
[HI, [49|. For these measurements we used the excited 
level l g as |e) and the "si" levels as \v) of the a 3 S+ 
state. The transition matrix element ci varies from about 
(0.2 to 33)MHz/Vw cm . In terms of a dipole moment, 
(er) = ^2(w)/v / ^2 x fry/eoc/2, this corresponds to (0.05 
to 8.0)xl0- 30 Cm. 

In comparison, the amplitude for the transition be- 
tween |i) and |e) is fii/V7i =0.4MHz/ VWcm" 2 . We 
determine Cl\ from resonant excitation, measuring how 
quickly |i) molecules are lost for a given laser intensity 
I\ . The data in Fig. [9] can be used to fix the position 
of the 1 3 E+ potential relative to the a 3 S^ potential ap- 
plying the Franck-Condon approximation. 



VI. CONCLUSION AND OUTLOOK 

In this paper we demonstrate high resolution dark state 
spectroscopy of ultracold 87 Rb2 molecules. We are able 
to resolve the vibrational, rotational, hyperfine and Zee- 
man structure of the lowest a 3 S+ potential with an ab- 
solute accuracy as high as 30 MHz. We find that the hy- 
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TABLE III: Scattering lengths (in Bohr radius, a = 0.5292 x 
10~ 10 m). ai owest gives the scattering length of the energeti- 
cally lowest hyperfine state at B z — 0. 



isotope 


& singlet 


^triplet 


^lowest R-b + Rb 
(/,"»/) + if,m f ) 


87/87 
85/85 
87/85 


90.35 
2720 
11.37 


99.04 
-386.9 
201.0 


100.36 (1,1) + (1,1) 
-460.1 (2, 2) + (2, 2) 
229.4 (1,1) + (2, 2) 



pcrfine structure is weakly dependent on the vibrational 
level, which also influences predictions of Feshbach reso- 
nances if one would like to reach accuracies on the order 
of 0.1 G. 

By optimizing mainly the triplet Born-Oppenheimer 
potential we obtain a model that can quite accurately 
predict all vibrational levels in the a 3 E + and X 1 S+ po- 
tentials for different ranges of rotational quantum num- 
bers. After reassi gnin g recent data from Fourier trans- 
form spectroscopy [17| we extend the range of applicabil- 
ity to rotational states as high as N — 70 for the a 3 S+ 
potential. Because of this reassignment also the molec- 
ular parameters like the equilibrium internuclear separa- 
tion R e or the dissociation energy D e have changed sig- 
nificantly compared to values reported in [l7|. The new 
values are given in Table |H] For rotational levels with 
N < 4 the model calculations for triplet levels of any 
vibrational level should have a precision similar to that 
of our measurements, i.e. about 30 MHz. For higher 
N this precision will increasingly degrade resulting from 
the reduced accuracy of the data from [ijj of about 300 
MHz for N = 70. Compared to [H the potential of 
the ground state X 1 ^ is significantly improved close 
to the atomic asymptote by including data on the mixed 
singlet-triplet levels of this study and data on Feshbach 
resonances from various other sources. We recommend 
the derived potential function for further use. It can 
predict the deeply bound levels with an accuracy (about 
50 MHz), comparable to that of the Fourier transform 
spectroscopy in 39] . The asymptotic levels are accurate 
on the order of a few MHz or better as their position 
is determined by the precisely measured location of Fes- 
hbach resonances. We improved the description of the 
large set of 87 Rb Feshbach resonances as compared to 
ref. [l9[ where the original calculation was based only on 
a few selected resonances and a derived asymptotic form 
of the two potentials. In [l9| deviations as high as 2 G 
appeared and for example the resonances of asymptote 
/a = l,mA = + /s = l,rag = 1 were all calculated 



systematically too high. In the present model these dis- 
crepancies disappear. The average deviation over all 46 
resonances, given in [i"9rl22l | . is about 0.15 G [5(|. If one 
removes the R-dependence of the hyperfine coupling as 
given in Eq.[8l deviations in the order of 0.5 G appear 
varying according to the different / levels which corre- 
late to the Feshbach resonances. To improve the overall 
fit, one probably needs to include many-body effects of 
the collision process for detecting the resonances. The 
overall improvement of the long range behavior also al- 
lows the whole set of scattering lengths for the three iso- 
topologues to be calculated. Table Hill gives a selection of 
these calculations. The error of these data will be in the 
last digit shown. These values are consistent with earlier 
publications but more precise and in general internally 
consistent between the isotopologues. 

Our model will be valuable for our planned collision 
experiments of ultracold molecules which have to be pre- 
pared in well defined quantum states. The singlet-triplet 
mixing that we observe is of interest for a proposed pre- 
cision experiment to measure the time dependence of the 
electron to proton mass ratio [HI (see also [j| [l(|). The 
precision of future measurements would be mainly lim- 
ited by the accuracy of the wavemeter and could be im- 
proved by one or two orders of magnitude. Such a high 
precision could be used, e.g. to test the fundamental lim- 
its of the coupled-channel model, where we assume the 
validity of (1) the Born-Oppenheimer approximation, (2) 
using Zeeman terms with atomic parameters only, (3) a 
limited functional dependence in R of the hyperfine and 
spin-spin interaction and (4) neglecting quadrupole hy- 
perfine coupling. 
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